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Multipole expansions depend on the coordinate system, so that coefficients of multipole moments 
can be set equal to zero by an appropriate choice of coordinates. Therefore, it is meaningless to 
say that a physical system has a nonvanishing quadrupole moment, say, without specifying which 
coordinate system is used. (Except if this moment is the lowest non-vanishing one.) This result 
is demonstrated for the case of two equal like electric charges. Specifically, an adapted coordinate 
system in which the potential is given by a monopole term only is explicitly found, the coefficients 
of all higher multipoles vanish identically. It is suggested that this result can be generalized to 
other potential problems, by making equal coordinate surfaces coincide with the potential problem's 
equipotential surfaces. 

I. INTRODUCTION 

Multipole expansions are a standard mathematical tool in both physics research and teaching (e.g., Ref.— ). For 
example, students of electromagnetism often calculate the electrostatic potential of a point electric charge when the 
coordinate system is not centered on the latter—^. In mathematical methods courses this problem often motivates, 
and serves as a means to the introduction of Legendre functions through its generating function 
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(The alternative approach to the introduction of the Legendre functions is through the angular equation obtained from 
the separation of variables of the Poisson equation in spherical symmetr y 4 ' 6 ' 11 ! 12 .) As the generating function is equal 
to the inverse distance between two points (in three-dimensional Euclidean space) — i.e., l/|x — x'| = g(r' /r, cos#), 
where x,x' are the two points — one can readily superpose solutions to obtain the field of an arbitrary distribution of 
point charges. Here, r, r' are their corresponding distances from the origin (so that r' < r) and 9 is the angle between 
the vectors x and x— . The electrostatic potential then is expressed as an infinite sum over Legendre functions (see, 
e.g.fi), with all multipoles present (both even and odd). 

The next problem in a typical mathematical methods course may often be to calculate the potential of a pair of 
charges, either of opposite signs ("electric dipole") or of the same sign. In either case, half of the coefficients in the 
series vanish, and one ends up with a series over all odd or even terms, respectively. Specifically, in the latter case 
(two equal charges q), the potential is 
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where r> (r< ) is the greater (smaller) of the evaluation point and the location of the charge from the center of the 
coordinate system. (We assume here, without loss of generality, that the charges are on the z-axis, so that the 
potential is azimuthal.) Notably, the potential includes a monopole term, a quadrupole term, a hexadecapole term, 
and so on. The different multipoles are used here describe the moments of a given charge distribution. Specifically, 
the quadrupole moment equals 2qr^ / r> . When r> 3> r< , the potential approximately equals the total charge divided 
by the distance to the evaluation point, which is the known result for a point particle at the origin of the coordinate 
system. Indeed, under this condition the separation of the charge from the origin is very small (compared with 
the distance to the evaluation point), so that the potential is indeed expected to be just such. However, when this 
condition is not satisfied, or indeed at finite distances, we find that the potential needs to be corrected, the correction 
terms being the multipole moments. 

Mathematical methods texts often discuss next how to arrange the charge distribution so that one, or more, of 
the multipole moments vanishes (in a given coordinate system)^. As the electrostatic problem is linear, one may 
superpose charges so that the lowest non- vanishing multipole moment may be as desired, e.g., one may obtain an 
electric dipole by having a charge q at z — a and a charge —q at z = —a thus canceling the monopole term, and the 
lowest non-vanishing moment is the dipole. An interesting point that is usually not highlighted is that the multipole 
expansion itself is coordinate dependent, so that for a given charge distribution one can choose adapted coordinated in 
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which certain multipole moments vanish. Specifically, we will show here in detail how, for the aforementioned example 
of two equal, like, static electric charges, one can find an adapted coordinate system in which all the multipoles above 
the monopole vanish. In such a coordinate system the potential is a function of one coordinate only, as the potential 
includes only a monopole term, which highlights the simplicity and elegance of the presentation. We see then that 
the complex multipole structure of the potential |T]) is a consequence of the choice of coordinates. We can find a 
coordinate system in which the potential equals the total charge divided by a coordinate that plays, in some loose 
sense, the role of "effective distance." It is only the lowest non-vanishing moment that cannot be nullified by a 
coordinate transformation. 

The price to pay for this simplicity and elegance is in the form of a non-trivial coordinate transformation from, 
say, cartesian coordinates to the adapted ones, and in topologically non-trivial coordinates. Arguably (and only 
semi-seriously), there is "conservation of complexity," which one cannot escape. However, several useful tools of 
mathematical physics are used in order to find this adapted coordinate system, which makes it an interesting problem, 
in addition to leading to deeper understanding of multipole expansions and their meaning. 



II. THE TWO-DISTANCE COORDINATES 



Let us first introduce a convenient coordinate system, from which we can transform to the desired adapted coor- 
dinates. Our strategy is to look for a coordinate system that is adapted to the potential in space. As the potential 
depends only on the distances to the two charges, a good starting point would be to base a coordinate system on 
the two distances from the two charges, so that equal coordinate surfaces would be spheres centered on either charge. 
The intersection of two such spheres spans a circle, and an azimuthal coordinate can then uniquely specify a point in 
space. To construct this coordinate system, position the charges along the x-axis, say, at x = ±a. Define the distance 
of the evaluation point of the field from the two centers: 



L := y/(x + a) 2 +y 2 + z 2 R := y/ (x - a) 2 + y 2 + z 2 . (2) 
These would be two coordinates, and the third one is defined by a rotation about the x-axis, namely, 



d> := tan- 1 (j j . (3) 

The coordinate transformation from the cartesian coordinates x, y, z to the coordinates L, R, $ is a regular coordinate 
transformation, as is evident from the explicit expressions for the inverse transformation: 
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The coordinates R, L, $ are the Two-Distance Coordinates, or 3D-TDC. We will find these coordinates to be a 
convenient starting point in our search for the adapted coordinate system, although they are unadapted to the 
potential problem themselves. The (spatial) metric in the 3D-TDC is given by 

= ( ^ +dR ^- ^l 2 ^^ ^ dR + ° 2 M2 ' (7) 

where 
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Q = R i6^ = 16^ =V +z . (8) 

The metric ([7]) is evidently non-diagonal, as it includes the "cross term" proportional to dL dR. Non-diagonal metrics 
are of much importance in physics and mathematical physics. Nevertheless, typical mathematical method courses 
do not discuss them. This metric is singular at the origin (half way between the two charges) — as can be seen from 
the vanishing of Q 2 there — but as this is flat Euclidean space, this singularity is a coordinate singularity. Figure [1] 
displays the 3D-TDC system. 



FIG. 1: The 3D-TDC: shown is the x — y plane. The two centers are at x = a,y — and x = — a.y = 0. The field evaluation 
point is at p. The coordinates L and R are defined by the distance of p from the two centers. The third coordinate <£> is 
obtained by rotation about the x-axis. 



III. SEARCHING FOR THE ADAPTED COORDINATES 



To search effectively for the adapted coordinate system, let us revisit regular spherical coordinates, and see why 
they are so effective in the description of the electric field of a single point charge. The main feature of the electric 
potential of a point charge is that equipotential surfaces are concentric spheres. The latter are the loci of all equidistant 
points from the center. These two statements imply that anywhere on an equipotential surface the value of one of 
the coordinates (the radial coordinate) is constant. To find the potential at some evaluation point, we only need to 
know on which equipotential surface it lies. For this reason, the potential is a function of only one coordinate, because 
only one coordiante is required to identify an equipotential surface. After the radial coordinate is chosen, one can 
readily find the other coordinates: the polar coordinate 9 is found by requiring that it is everywhere orthogonal to 
constant-r surfaces, and the azimuthal coordinate is obtained by rotation of the r — 9 plane. Spherical coordinates 
are therefore adapted to the potential problem of a single charge at their center, and the potential is described in full 
by a monopole term only. We show this derivation in detail in Appendix A. 

In searching for the adapted coordinates we follow a similar path: First, define the coordinate x so that equal- x 
surfaces are also equipotential surfaces of the static 3D Coulomb problem. Because of the problem's linearity, we 
may superpose solutions of two single-charge problems, for which the potentials (of unit charges) are 1/L and l/R, 
respectively, i.e., the total potential is 1/L + l/R. In fact, any monotonic function of 1/L + l/R is a good choice for 
this coordinate, i.e., \ = x{^/L + l/R). Specifically, we choose \ so that at great distances (L, R 3> a), the coordinate 
X behaves like the radial coordinate r, because in this limit the potential is almost that of a point particle (of charge 
2q). That is, we choose 
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Notice that the coordinate x is non-trivial: it undergoes a topology change at the critical surface x = a - For x < a 
equal- x surfaces are doubly connected, and for x > a they are singly connected. Figure [2] shows the equal-x surfaces. 

To find the coordinate we require that it is everywhere orthogonal to equal-x surfaces (which are also the 
equipotential surfaces). Namely, equal-0 lines will be the electric field lines. The condition that two curves are 
orthogonal is that their gradients are. The reason why we need to require the condition of the gradients is that inner 
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FIG. 2: The adapted coordinate system in the x — y plane. Rotation about the x-axis yields the full 3D system. The closed 
contours are the equal-x surfaces. The "figure eight" curve is the critical surface across which the connectedness of the x 
coordinate undergoes a topology change. The open curves are the equal-O lines. 



product is a map from two vectors to a scalar. Orthogonality of two vectors is defined by the requirement that their 
inner product vanishes. The way to create a vector from a scalar is by calculating its gradient. Specifically, satisfies 

Vx • V9 = 0, (11) 

which in Cartesian coordinates is 

dx dQ_ + dxdQ + d\ dQ_ _ Q ^ 
dx dx dy dy dz dz 

The most straighforward way to try and solve the partial differential equation for O (recall that \ is known; the 
unknown in Eq. (fTTj) is 6) is to try and solve Eq. (jTSJ) in cartesian coordinates. This turns out to be a highly non- 
trivial thing to do. In Appendix B we discuss this equation. The easier way to solve Eq. (jlip is by realizing that it is 
a scalar equation, so that we can express the terms on either hand side in whatever coordinates we find convenient, 
and the equation will retain its form. Equation pip becomes particularly simple when expressed in 3D-TDC. The 
reason why this could be expected is that space is symmetrical under the exchange of the two centers. It is therefore 
reasonable to expect that geometrical expressions would be simpler when expressed in terms of the distances to the 
two centers (and not other points such as the origin of the cartesian grid), the only "natural" quantities that the 
metric structure of space depends on. 

We can express Eq. (fTTj) in 3D-TDC in two way: first, we can transform Eq. (fT!?j) or we can use the metric (J7J) to 
write g lJ x,i®,j = 0, a comma denoting partial derivative. Either way, we find that Eq. (|11[) transforms to 
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whose solution is 

f R - T, . „ 1 

(14) 



6(L )J R) = e{^ [(i? + L) 2 -4a 2 ]} 



This solution is easy to verify, or even derive using software such as Maple or Mathematica. Notice that the solution 
is any function of a certain combination of L, R, with a acting as a parameter (the only length scale in this problem). 
Indeed, any function of this argument is everywhere orthogonal to equal- % surfaces. To make our solution useful, we 
next need to choose judiciously which function. Specifically, we make our choice by requiring that the asymptotic 
properties of the solution for O will coincide with the regular polar coordinate 9. This way, our coordinate system 
will asymptotically approach regular spherical coordinates as it should, as at great distances the separation of the two 
centers becomes negligible. 
At very great distances, 

^ [(i* + Lf - 4a 2 ] -> -8a- + 0(y- 3 ) . 
RL y 

We next require that at that limit a coordinate coincides with the regular polar coordinate 9 := tan _1 (y/x), which 
motivates us to choose 



B(L,R) = tan- 
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This guarantees that at very great distances behaves similarly to the regular polar coordinate 9. However, the 
coordinate is still not quite what we need, because it is discontinuous: Consider a semicircle at constant large 
distance from the center of the cartesian coordinate system (the center is half way between the two charges) , starting 
on the positive x-axis and going through the upper half-plane to the negative x-axis. The value of will vary from 
7r/4 to 7r/2, jump discontinuously to — tt/2 crossing the y-axis, and then change to — 7r/4. Define then 

0(L, R) := | - TT sgn(L - R) + 2Q(L, R) . (16) 

The coordinate is continuous in the upper half-plane, and it is everywhere orthogonal to x- It's range is from 
to n. Notice that sgn(L — R) = sgn(x) = sgn(7r/2 — 0), and that we only need to define in the upper half-plane, 
because we rotate the coordinates about the x-axis. 

We now have the coordinate system 0,$), which is the Two-Center Bi-Spherical coordinates system (TCBS). 
The metric in these coordinates is given by 

, 2 = 1 (L + Rf 2 

° A {L + Rf- 3LR(L + R) 2 - ALRa 2 X 

, 1 WL-W-HL' + V-WrWHL-mL+msi. 



16 • 256 a 4 Q 2 [(L + R) A - 3LR(L + R) 2 - ALRa 2 ] 

+ Q 2 d$ 2 , (17) 

where L, R are implicit functions of Xi 0- The metric ([T7|) is manifestly diagonal. Given %, 0, one needs to solve the 
following quintic to find R(x, 0) explicitly: 

R 5 ~ X R 4 + Mh^X - a) R 3 - Aaxih^x - 2a) R 2 + ax 2 {hr x x -5a)R + a 2 X 3 = (18) 

where 



h := tan = tan — 
2 



0-| + 7rsgn(| -0 



As is well known from Abel's Impossibility Theorem 1 , it is impossible to solve a general quintic equation in terms of 
radicals. However, solutions in terms of hypergeometric functions^ or Jacobi Theta functions^ are always possible. 
Numerical solutions are of course easy to find (e.g., using the Newton-Raphson method). Notice, that the Fundamental 
Theorem of Algebra guarantees a real solution. It is immediately clear that the coordinate system is singular at the 
two centers and at the origin of the cartesian coordinate system. For example, Q 2 vanishes at all three singular 
points, so that gee diverges, and vanishes. As space is 3D flat Euclidean space, we know that this singularity is a 
coordinate singularity, and not a genuine geometrical one. Notably, the Jacobian of the transformation from cartesian 
to 3D-TCBS coordinates is regular. 

Given the solution for R(x, 0) we can readily find L(x, 0): 
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IV. THE POTENTIAL PROBLEM IN THE ADAPTED COORDINATES 

To express the potential problem in the newly found TCBS coordinates we first express the Laplacian in these 
coordinates. In any coordinate system the Laplacian of a scalar field \& is given by 

v 2 * = ( Vx • v%) * xx + ( ve • ve) *, e + ( v$ ■ v$) 

+ 2( V% • V6) *, x e + 2( Vx • V4>) # x# + 2( V9 • V$) #, e * 

+ (V 2 x)*, x + (V 2 8)*,e + (V 2 $)*, <[ ,. (20) 

The coefficients in this expression can readily be calculated explicitly from the above expressions. Specifically, we find 

L 4 + R 4 + LR(L 2 + R 2 - 4a 2 ^ 
LR{L + Rf 



V 2 X^4 - ^ -3 (21) 



^2^ 128 a 3 {L-R)G 

V ft — : i (221 

[16(L + i?) 2 a 4 + 8(L 4 -10L 2 i? 2 + i? 4 )a 2 -( J L-i?) 2 ( J L + i?) 4 ] 2 1 ' 

where 

G = 64a 6 - 16(5L 2 + 2LR + 5i? 2 )a 4 + 4(7L 4 + 28L 3 i? + 26L 2 i? 2 + 28LR 3 + 7R 4 )a 2 

- (L - R) 2 (3L 4 + 20L 3 R - UL 2 R 2 + 20LR 3 + 3R 4 ) , (23) 

V 2 $ = (24) 

Vx • V X = 9 XX (25) 

V6 • V6 = g ee (26) 

V$ • V$ = <?** (27) 

Vx • Ve = (28) 

Vx • V$ = (29) 

ve • v$ = , (30) 

the last three relations resulting from the TCBS coordinates being orthogonal. Notice that $ is harmonic. 

We are now in a position to show that in the 3D-TCBS the solution for the potential problem of two equal like 
charges is given by a monopole term only. Specifically, we show that 

V 2 (£)=0 (x^O). (31) 

This relation can be directly verified by substitution in Eq. (|2"0)) : Clearly all the derivatives with respect to either Q 
or $ vanish, so that Eq. (|20| reduces to 

v*(I)= 9 «(I) 

\ A / \ A / , xx \ A 

2 V 2 y 

X 3 X 

after direct substitution of g xx from Eq. (jTTJ) and using Eq. pTj) . Equation (I3L[) is analogous to V 2 (l/r) = (except 
at r = 0) in regular spherical coordinates. At x = (x = ±a, y = z = 0) the Laplacian of x no longer vanishes, so that 
the global problem is described by Poisson's equation, V 2 ^ = SnqJ^ 1 S(x), where J is the Jacobian determinant 
for the transformation from cartesian to 3D-TCBS coordinates. Comparing Poisson's equation with Eq. (|3"Tj) , we find 
the solution for the potential problem to be 

* = ^ , (33) 
X 

which is the desired solution. 



7 



V. 



CONCLUSIONS 



We showed, for the specific example of two equal like electric charges, how to find a coordinate system in which the 
electric potential is described by one coordinate, and for which it is given by a monopole term only. The derivation 
involved a large number of topics that are covered in mathematical methods courses, such as non-diagonal coordinate 
systems, coordinate singularities, quintic equations, multipole expansions, coordinate transformations, and potential 
theory. As such, it may serve as an instructive problem for such courses. In particular, it may be used to demonstrate 
the deep meaning of a multipole expansion, and its dependence on the choice of coordinates. Specifically, a multipole 
expansion can be made simple when the coordinates used are adapted to the potential problem. The reason why 
an infinite number of (even) multipoles are needed to describe the potential of two equal like charges using regular 
spherical coordinates JT]) is that the equal coordinate surfaces are very different from the equipotential surfaces. By 
choosing the two surfaces to coincide, we are able to eliminate all the higher moltipoles from the potential, and solve 
the potential problem using only the monopole term. The price to pay is in a form of more involved mathematics, 
which is however still in the range of the usual preparation of the usual Physics education programs. 



The author wishes to thank Richard Price and Anthony Hester for discussions, and Ross Cortez for checking some 
of the calculations. This work was supported in part by a minigrant from the UAH Office of the Vice President for 
Research. 



In this Appendix we apply our method to a point charge, and derive the adapted coordinates, which are just the 
regular spherical coordinates. This derivation may serve as a pedagogic illustration of the method, applied to a trivial 
situation for which the solution is well known. 

As discussed above, for the case of a single point particle, the radial coordinate r describes equipotential surfaces. 
We therefore choose the radial coordinate r as one of the adapted coordinates, or \s '■= T. To find the second 
coordinate, we require that it is orthogonal to r. Therefore, it satisfies Eq. (11), i.e., Vx s ■ V0 S = 0, or, in 
cartesian coordinates in the x-y plane, (dr/ dx)( d@ s / dx) + ( dr dy)/( d<d s dy) = 0. Recalling that drj dx = x/r 
and drj dy — y/r, this equation becomes x ( d<d s / dx) + y ( dQ s / dy) — (except for the origin which is a coordinate 
singularity), whose solution is O s = Q s (y/x). Next, we require that at large distances this coordinate coincides with 
the regular azimuthal coordinate 8 := taxi (y/x), so that we find that everywhere O s := 9. 

Rotating the coordinates by an angle <f> about the cc-axis yields x' = x = r cos 8, y' = y cos0 = r sin# cos</>, and 
z' = y sin</) = r sin# sin</>. The regular spherical coordinates are obtained by renaming the cartesian axes x" = y' , 
y" = z', and z" = x': x — r sin# cos</>, y = r sin (9 sin^, and z = r cos 8, after dropping the primes for conventionality. 



Equation (|12p can be readily written explicitly. In the x — y plane, the unknown is 0(x, y), so that Eq. (|12[) becomes 



This partial differential equation is non-trivial to solve symbolically. Standard techniques, e.g., separation of variables, 
prove to be ineffective. Even PDE solvers such as Maple and Mathematica, using up to 8 gigabytes RAM for a couple 
of weeks on a PowerMac G5 were unsuccessful in solving this equation. Indeed, having found the solution (|I4p using 
the 3D-TDC it is evident why straightforward attempts to solve this equation failed. It should be noted, however, 
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APPENDIX A: "DERIVATION" OF SPHERICAL COORDINATES 



APPENDIX B: TRYING TO SOLVE EQ. (TT2|) DIRECTLY 
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that numerical solutions for this equation can readily be found, if appropriate boundary conditions are specified. 
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